Search for: All records

Context.One of the most important open questions in planet formation is how dust grains in a protoplanetary disk manage to overcome growth barriers and form the ∼100 km planet building blocks that we call planetesimals. There appears to be a gap between the largest grains that can be produce by coagulation, and the smallest grains that are needed for the streaming instability (SI) to form planetesimals. Aims.Here we explore the novel hypothesis that dust coagulation and the SI work in tandem; in other words, they form a feedback loop where each one boosts the action of the other to bridge the gap between dust grains and planetesimals. Methods.We developed a semi-analytical model of dust concentration due to the SI, and an analytic model of how the SI affects the fragmentation and radial drift barriers. We then combined them to model our proposed feedback loop. Results.In the fragmentation-limited regime, we find a powerful synergy between the SI and dust growth that drastically increases both grain sizes and densities. We find that a midplane dust-to-gas ratio ofϵ ≥ 0.3 is a sufficient condition for the feedback loop to reach the planetesimal-forming region for turbulence values 10⁻⁴ ≤ α ≤ 10⁻³and grain sizes 0.01 ≤ St ≤ 0.1. In contrast, the drift-limited regime only shows grain growth without significant dust accumulation. In other words, planetesimal formation remains challenging in the drift-dominated regime and dust traps may be required to allow planet formation at wide orbital distances. 

Abstract The streaming instability (SI) is a leading mechanism for concentrating solid particles into regions dense enough to form planetesimals. Its efficiency in clumping particles depends primarily on the dimensionless stopping time (τ_s, a proxy for particle size) and dust-to-gas surface density ratio (Z). Previous simulations identified a criticalZ(Z_crit) above which strong clumping occurs, where particle densities exceed the Hill density (thus satisfying a condition for gravitational collapse), over a wide range ofτ_s. These works found that, forτ_s≤ 0.01,Z_critwas above the interstellar medium value (∼0.01). In this work, we reexamine the clumping threshold using 2D axisymmetric, stratified simulations at high resolution and with relatively large (compared to many previous simulations) domain sizes. Our main results are as follows: First, whenτ_s = 0.01, strong clumping occurs even atZ ≲ 0.01, lower thanZ_critfound in all previous studies. Consequently, we revise a previously published fit to theZ_critcurve to account for this updatedZ_crit. Second, higher resolution results in a thicker dust layer, which may result from other instabilities manifesting, such as the vertically shearing SI. Third, despite this thicker layer, higher resolution can lead to strong clumping even with a lower midplane dust-to-gas density ratios (which results from the thicker particle layer) so long asZ ≳ Z_crit. Our results demonstrate the efficiency of the SI in clumping small particles atZ ∼ 0.01, which is a significant refinement of the conditions for planetesimal formation by the SI. 

Abstract Ionization drives important chemical and dynamical processes within protoplanetary disks, including the formation of organics and water in the cold midplane and the transportation of material via accretion and magnetohydrodynamic flows. Understanding these ionization-driven processes is crucial for understanding disk evolution and planet formation. We use new and archival Atacama Large Millimeter/submillimeter Array observations of HCO⁺, H¹³CO⁺, and N₂H⁺to produce the first forward-modeled 2D ionization constraints for the DM Tau protoplanetary disk. We include ionization from multiple sources and explore the disk chemistry under a range of ionizing conditions. Abundances from our 2D chemical models are postprocessed using non-LTE radiative transfer, visibility sampling, and imaging, and are compared directly to the observed radial emission profiles. The observations are best fit by a modestly reduced cosmic-ray ionization rate (ζ_CR∼10⁻¹⁸s⁻¹) and a hard X-ray spectrum (hardness ratio = 0.3), which we associate with stellar flaring conditions. Our best-fit model underproduces emission in the inner disk, suggesting that there may be an additional mechanism enhancing ionization in DM Tau’s inner disk. Overall, our findings highlight the complexity of ionization in protoplanetary disks and the need for high-resolution multiline studies. 

Abstract The streaming instability (SI) is a leading candidate for planetesimal formation, which can concentrate solids through two-way aerodynamic interactions with the gas. The resulting concentrations can become sufficiently dense to collapse under particle self-gravity, forming planetesimals. Previous studies have carried out large parameter surveys to establish the critical particle to gas surface density ratio (Z), above which SI-induced concentration triggers planetesimal formation. The thresholdZdepends on the dimensionless stopping time (τ_s, a proxy for dust size). However, these studies neglected both particle self-gravity and external turbulence. Here, we perform 3D stratified shearing box simulations with both particle self-gravity and turbulent forcing, which we characterize via a turbulent diffusion parameter,α_D. We find that forced turbulence, at amplitudes plausibly present in some protoplanetary disks, can increase the thresholdZby up to an order of magnitude. For example, forτ_s= 0.01, planetesimal formation occurs whenZ≳ 0.06, ≳0.1, and ≳0.2 atα_D= 10⁻⁴, 10^−3.5, and 10⁻³, respectively. We provide a single fit to the criticalZrequired for the SI to work as a function ofα_Dandτ_s(although limited to the rangeτ_s= 0.01–0.1). Our simulations also show that planetesimal formation requires a mid-plane particle-to-gas density ratio that exceeds unity, with the critical value being largely insensitive toα_D. Finally, we provide an estimation of particle scale height that accounts for both particle feedback and external turbulence. 

Abstract Kuiper Belt objects (KBOs) show an unexpected trend, whereby large bodies have increasingly higher densities, up to five times greater than their smaller counterparts. Current explanations for this trend assume formation at constant composition, with the increasing density resulting from gravitational compaction. However, this scenario poses a timing problem to avoid early melting by decay of²⁶Al. We aim to explain the density trend in the context of streaming instability and pebble accretion. Small pebbles experience lofting into the atmosphere of the disk, being exposed to UV and partially losing their ice via desorption. Conversely, larger pebbles are shielded and remain icier. We use a shearing box model including gas and solids, the latter split into ices and silicate pebbles. Self-gravity is included, allowing dense clumps to collapse into planetesimals. We find that the streaming instability leads to the formation of mostly icy planetesimals, albeit with an unexpected trend that the lighter ones are more silicate-rich than the heavier ones. We feed the resulting planetesimals into a pebble accretion integrator with a continuous size distribution, finding that they undergo drastic changes in composition as they preferentially accrete silicate pebbles. The density and masses of large KBOs are best reproduced if they form between 15 and 22 au. Our solution avoids the timing problem because the first planetesimals are primarily icy and²⁶Al is mostly incorporated in the slow phase of silicate pebble accretion. Our results lend further credibility to the streaming instability and pebble accretion as formation and growth mechanisms.